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ABSTRACT 

The hyper luminous X-ray source HLX-1 in the galaxy ESO 243-49, currently the best in- 
termediate mass black hole candidate, displays spectral transitions similar to those observed in 
Galactic black hole binaries, but with a luminosity 100-1000 times higher. We investigated the 
X-ray properties of this unique source fitting multi-epoch data collected by Swift, XMM-Newton 
& Chandra with a disk model computing spectra for a wide range of sub- and super-Eddington 
accretion rates assuming a non-spinning black hole and a face-on disk (i = 0°). Under these 
assumptions we find that the black hole in HLX-1 is in the intermediate mass range (~2x 10 4 
Mq) and the accretion flow is in the sub-Eddington regime. The disk radiation efficiency is 
rj = 0.11 ± 0.03. We also show that the source does follow the Lx oc T A relation for our mass 
estimate. At the outburst peaks, the source radiates near the Eddington limit. The accretion rate 
then stays constant around 4 x 10~ 4 M Q yr _1 for several days and then decreases exponentially. 
Such "plateaus" in the accretion rate could be evidence that enhanced mass transfer rate is the 
driving outburst mechanism in HLX-1. We also report on the new outburst observed in August 
2011 by the Swift-X-ray Telescope. The time of this new outburst further strengthens the ~ 1 
year recurrence timescale. 

Subject headings: Galaxies: individual (ESO 243-49); Physical data and processes: Accretion, accretion 
disks, Black hole physics; Methods: data analysis; X-rays: individual (HLX-1) 

1. Introduction 
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Ultra Luminous X-ray sources (ULXs) are de- 
fined as off-nucleus extragalactic sources showing 
X-ray luminosity exceeding 3 x 10 39 erg s _1 assum- 
ing isotropic emission (see Roberts et al. 2007). 
Even if their nature is still in dispute, it is likely 
that their huge luminosity is produced by accre- 
tion of matter onto a black hole (BH). Three ex- 
planations for their nature have been considered. 
(1) ULXs may be X-ray stellar-mass BH binaries 
(BHBs) similar to those observed in our Galaxy, 
but in a more extreme version of the very high 
state (e.g. Remillard & McClintock 2006), the ul- 
traluminous or wind-dominated state (e.g. Glad- 
stone et al. 2009). In such states, the source would 
be able to radiate above the Eddington limit. It is 
still to be understood why this ultraluminous state 
is so rarely observed in X-ray binaries. (2) King et 
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al. (2001) proposed that the emission of ULXs is 
highly anisotropic. In this case, the requirement to 
have super-Eddington emission is alleviated. The 
nature of this anisotropy could be due to either ge- 
ometrically thick disks funnelling the X-ray pho- 
tons produced in the inner parts of the accretion 
disks (King 2009) or due to relativistic beaming 
of a jet or strong outflow. However, the discovery 
of several optical and radio nebulae around ULXs 
(e.g. Pakull & Grise 2008) and QPOs in M82 X- 
1 (Strohmayer & Mushotzky 2003) argue against 
strong beaming. (3) The final and the most ex- 
citing explanation is that some ULXs are accret- 
ing intermediate mass BHs (IMBHs) with masses 
ranging from ~ 100 M Q to ~ 10 5 M (e.g. Col- 
bert & Mushotzky 1999). The existence of such 
IMBHs will naturally alleviate the need for super- 
Eddington emission. It is clear now that the IMBH 
interpretation is not valid to explain the ULX pop- 
ulation as a whole (see e.g. Roberts 2007). How- 
ever, the most luminous ULXs, the so-called Hy- 
per Luminous X-ray sources (HLX - Gao et al. 
2003) with X-ray luminosities above 10 41 erg s _1 , 
are good candidates. Finding convincing evidence 
for the existence of IMBHs is important for the 
growth of supermassive BHs via mergers or accre- 
tion episodes (e.g. Micic et al. 2007), dark matter 
studies (e.g. Fornasa & Bertone 2008, but see also 
Bringmann et al. 2009), cosmology (e.g. Trenti 
& Stiavelli 2007) and gravitational wave detection 
(e.g. Matsubayashi et al. 2004, Amaro-Seoane & 
Santamara 2010). 

Farrell et al. (2009) reported the serendipitous 
discovery of a ULX candidate 2XMM J011028.1- 
460421, referred to hereafter as HLX-1, located in 
the outskirts of the edge-on spiral galaxy ESO 243- 
49 at a redshift of 0.0224 (Wiersema et al. 2010). 
From its maximum luminosity reaching ~ 1.3 x 
10 42 erg s _1 at peak and assuming that the source 
luminosity reached up to ten times the Eddington 
limit (L Edd ), Farrell et al. (2009) derived a BH 
mass of more than 500 M©. HLX-1 is so far the 
best candidate to harbour an IMBH. Apart from 
its extreme luminosity which has been observed 
many times over the past 3 years by different X- 
ray satellites (Swift, XMM-Newton & Chandra), 
HLX-1 is unique amongst other ULXs because it 
is the only one for which clear spectral hysteresis 
similar to those observed in Galactic BHBs (GB- 
HBs) are seen. Indeed, Servillat et al. (2011) using 



XMM-Newton, Chandra and Swift data confirmed 
the results presented in Godet et al. (2009) that 
HLX-1 underwent transitions from the high/soft 
state to the low/hard state. The Swift-X-ray Tele- 
scope (XRT) lightcurve covers so far four out- 
bursts from 2008 to 2011 (hereafter P0, PI, P2 
and P3 - see Fig. [T]). While we had only one ob- 
servation in October 2008 during the P0 outburst, 
the PI and P2 outbursts display well-sampled 
Fast-Rise-Exponential-Decay (FRED)-like tempo- 
ral profiles. Thanks to our dedicated Swift-XKT 
ToO, we caught the rise of the P3 outburst from 
15 th August 2011 (Godet et al. 2011). A precur- 
sor could be seen prior to the P3 outburst peak 
in Fig. [TJ Two reflare events are also visible dur- 
ing the PI and P3 outbursts. The outbursts ap- 
pear to be separated by a recurrence timescale of 
nearly a year. Recently, Lasota et al. (2011) in- 
terpreted the X-ray light-curve as the result of en- 
hanced mass-transfer rate onto a pre-existing ac- 
cretion disk around an IMBH, when an asymptotic 
giant branch star orbiting along an eccentric orbit 
with a period of ~ 1 year is tidally stripped near 
periastron. 

In order to further investigate the nature of 
this unique source, it is essential to put some con- 
straints on the accretion flow and the BH mass. 
Lower and upper limits from radio, optical and X- 
ray observations have been derived (see Wiersema 
et al. 2011, Servillat et al. 2011 and Webb et al. 
2011). Davis et al. (2011) using their relativistic 
accretion disk model bhspec to fit X-ray spectra 
when the source was in various spectral states put 
some constraints on the BH mass. They showed 
that the inclination (i) has a strong influence on 
both the mass and the spin (a*) of the BH de- 
rived by the model. Due to degeneracies in their 
best-fit parameters, they were only able to give a 
range of possible BH masses from 3000 M© (where 
the limit corresponds to i = 0°, T Lx = 0.7 and 

^Edd 

a* = -1) to 3 x 10 5 M© (where the limit cor- 
responds to i — 90° and a* = 0.99). However, 
the inclination is likely to be less than 60 — 70° 
due to the lack of observed eclipse dips assum- 
ing HLX-1 is a binary system. These results were 
obtained neglecting the effects of radial advection 
that are important at luminosities above the Ed- 
dington limit L Edd with L Edd = 1.3 x 10 38 

erg s . Radial advection is a mechanism enabled 
to stabilize disks in the super-Eddington accretion 
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regime (see e.g. Abramowicz et al. 1988). This 
regime is reached when the accretion rate (M) is 
larger than MEdd = r/ 1 L ^z d w itti c the speed of 
light and n the radiation efficiency. In this regime, 
the vertical disk structure could be geometrically 
thick (H/R > 0.5 with H and R the scale height 
and the radius of the disk, respectively). How- 
ever, for a range of moderate super-Eddington M- 
values there are still disk solutions for which the 
vertical disk structure can be considered as rela- 
tively thin. They are often referred to as slim disk 
solutions (Abramowicz et al. 1988). For large ac- 
cretion rates (rh > 10 with rh — ¥- c ), electron 
scattering (opacity and comptonization) have sev- 
eral effects on the emergent spectra: i) they can be 
highly distorted and no longer look like a multi- 
color black-body (BB) spectrum (e.g. Kawaguchi 
2003, hereafter K03); ii) the L oc T 4 relation with 
T the colour temperature at the inner disk radius 
for standard a-disk (Shakura & Sunayev 1973) is 
no longer valid (see Fukue 2000, K03). In addi- 
tion, slim disks may extend beyond the innermost 
stable circular orbit (ISCO) due to non-negligible 
pressure support resulting in higher disk tempera- 
tures (e.g. Abramowicz et al. 1988, 2010; Watarai 
et al. 2000). 

In this paper, we make use of the K03 disk 
model in order to: i) investigate the accretion 
flow and BH properties in HLX-1 through spec- 
tral fitting of a multi-epoch and multi-instrument 
dataset using a model that computes a wide range 
of sub- and super-Eddington accretion disk solu- 
tions; ii) investigate whether radial advection and 
electron scattering have important effects on the 
emergent spectra; hi) compare our results with 
those performed on less luminous ULXs using the 
same model that always favoured critical or super- 
Eddington accretion onto a stellar mass BH (e.g. 
Yoshida et al. 2010, Foschini et al. 2006, Vier- 
dayanti et al. 2006, Okajima et al. 2006). 

The paper is organised as follows. In Section 2, 
we describe the data used to perform the spectral 
analysis, as well as the data reduction. Section 3 
presents some possible observational evidence for 
effects of radial advection and/or electron scatter- 
ing as well as the detection of possible X-ray lines. 
In Section 4, we present the fitting results obtained 
with the K03 disk model. In Section 5, we discuss 
the implications of the fitting results on the BH 
mass, the accretion rate and the disk structure. 



The main conclusions are given in Section 6. 

2. Observation log & Data reduction 

We consider X-ray data coming from 3 observa- 
tories: XMM-Newton, Chandra (ACIS) and Swift. 
The data used in this study cover the different 
spectral states observed during the outbursts of 
HLX-1 (see Table EJ). We used the same XMM- 
Newton and Chandra spectra as those presented 
in Servillat et al. (2011). The nomenclature used 
here is the same as that in Servillat et al. (2011) 
and Farrell et al. (2009) for the XMM-Newton 
data (see Table [3]). For the XMM-Newton data, 
we only used the EPIC-pn data because it offers 
the best statistics for a given observation. Please 
refer to Servillat et al. (2011) for the details of the 
data reduction on the XMM-Newton and Chan- 
dra data. In order to investigate for the presence 
of X-ray lines during the outbursts, we also used 
the Reflection Grating Spectrometer (RGS) data 
from the XMM2 observation (see Table 13]). The 
RGS data have been reduced using the standard 
XMM-Newton Science Analysis Software (SAS) 
task rgsproc and the most recent calibration files 
released in February 2011. For the spectral analy- 
sis, we only considered the first order of the RGS 
data. 

All the Swift-XBT Photon Counting data were 
processed using the HEASOFT v6.11 and the lat- 
test calibration files (CALDB version 3.8). This 
new CALDB includes a new gain file enabled to 
correct the data for charge traps that accumu- 
late on the CCD due to radiation damage (see 
Pagani et al. 2011). These traps can induce 
some energy offsets. The data were processed us- 
ing the tool xrtpipeline vO.12.6. We inspected 
all segments of data to search for epochs when 
the background was enhanced by the presence 
of hot pixels induced by a high CCD tempera- 
ture (T > —55°). For these time intervals, we 
checked whether there were some hot pixels in 
the extraction regions used to extract the spec- 
tra. Thus, we excluded the data from segments 
00031287125 (0.8 ks) and 00031287129 (0.2 ks), be- 
cause the data were severely contaminated by hot 
pixels due to a high (T > -51°C) CCD tempera- 
ture. We also excluded the first orbit (~ 0.6 ks) of 
segment 00031287023 for the same reason. These 
latter bad data induce an artificially high count 
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rate in the X-ray light-curve at MJD = 55209 
(during the PI outburst - see Fig. [T]). However, 
the re-flare event seen around that time is not due 
to hot pixel contamination. We used the grade 
0-12 events, giving slightly higher effective area 
at higher energies than the grade events, and 
a 20 pixel (47.2 arcsec) radius circle to extract 
the source and background spectra using xselect 
v2.4b. The background extraction region was cho- 
sen in order to be close to the source extraction re- 
gion and not to contain any XMM- Newton sources. 
The ancillary response files were generated us- 
ing xrtmkarf vO.5.9 and exposure maps gener- 
ated by xrtexpomap vO.2.7. The response file 
swxpcotoi2S6_2ooioioivoi2.rmf is used to fit the 
spectra. To extract the Swift-XRT spectra, we di- 
vided the data over different intensity ranges for 
the outbursts P0, Pf and P2. The work on the 
first three outbursts being completed before the 
most recent outburst P3, we decided to compare 
the latest outburst with the previous three by di- 
viding the data over time instead. Table [4] sum- 
marises the nomenclature used for the Swift-XRT 
data. 

3. Spectral analysis 

All the spectra were grouped at a minimum of 
20 counts per bin to provide sufficient statistics to 
use the \ 2 minimisation technique. The XMM3, 
Swift-XRT and Chandra spectra were fitted in the 
0.3-10 keV energy range within xspec vl2.7.0 (Ar- 
naud 1996), while the XMM1 and XMM2 spectra 
were fitted in the 0.2-10 keV energy range. The 
Chandra data being moderately piled-up were fit- 
ted using the pileup model (Davis 2001) within 
xspec with a frame time fixed at 0.8 s (see Servil- 
lat et al. 2011 for more details). For each model, 
we determined the appropriate value of the grade 
morphing parameter. For the data used in this 
paper, the statistical errors are dominant over the 
instrumental systematics. The RGS1 and RGS2 
spectra were binned at twice the resolution of the 
instrument (AA = 0.2 A). Since with this choice 
of binning there are fewer than 20 counts per res- 
olution bin, we used the C-statistic (Cash 1979) 
available within xspec for the spectral fit. 

The total column density was fixed at the best- 
constrained value of Nu = 4 x 10 20 cm -2 from the 
XMM2 observation (see Farrell et al. 2009). The 



Galactic absorption column in the direction of the 
source is equal to 1.8 x 10 20 cm -2 (e.g. Kalberla et 
al. 2005). For all fits, the absorption is modeled 
using the wabs photo-electric absorption model. 
The source redshift being 0.0224 (Wiersema et al. 
2010), we adopted a source distance of dj, = 95 
Mpc using the cosmological parameters from the 
WMAP5 results (H = 71 km s" 1 Mpc" 1 , Cl M = 
0.27 & =0.73). All the errors quoted below are 
given at a 90% confidence level for one parameter 
of interest (i.e. A X 2 = 2.706). 

3.1. Hardness-Intensity diagram 

We fitted the Swift-XRT spectrum using a 
multi-color disk diskbb and/or powerlaw model. 
This phenomelogical model is often used to fit the 
spectra of GBHBs. In that context, the power- 
law component is often interpreted as emission 
produced by a corona of hot electrons located in 
the inner regions of the accretion disk, while the 
diskbb component is interpreted as the emission 
from the accretion disk (see Remillard & McClin- 
tok 2006). The photon index of the powerlaw 
was tied together for the Scm,i, 2,3,4,5 spectra, 
because we found T-values consistent within the 
errors when fitting the spectra individually. The 
powerlaw normalisation and the parameters of the 
diskbb model were left free to vary independently 
between spectra. The 5*41,2 spectra are well fitted 
by an absorbed diskbb model. We did not obtain 
a good fit (x 2 /dof = 14.5/8) using an absorbed 
diskbb model for the St3 spectrum. There are 
residuals left at high energies. When adding a 
powerlaw component with F free, the T parame- 
ter was not constrained. So, we decided to freeze 
this parameter at the same value as found for the 
Scro,i,2, 3,4,5 spectra (r = 2.2). The fit is im- 
proved (A% 2 = 5 for 1 dof i.e. an improvement 
at a > 2 a significance level). For each spectrum, 
we then computed the hardness ratio (HR) of the 
0.3-1.5 keV observed flux over the 1.5-10 keV flux. 
This is different from the Swift hardness ratio pre- 
sented in Godet et al. (2009) and Servillat et al. 
(2011) that was based on count rate. The HR 
values, the best-fit spectral parameters and the 
0.2-10 keV unabsorbed luminosity (L = Air d\ F 
with F, the 0.2-10 keV unabsorbed flux) for each 
spectrum are summarised in Table [5j 

Fig. [2] shows the Swift-XRT hardness ratio ver- 
sus the 0.2-10 keV unabsorbed luminosity. The 
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points from the Chandra (kT = 0.21 ± 0.01 keV, 
N = 37t"), XMM2 (kT = 0.17 ± 0.01 keV, 
N = 40±H, r = 2.1 ± 0.4, N r = l$t%i x 10" 5 
ph cm- 2 s- 1 keV" 1 ) and XMM30 (r = 2.2t°; 4 
and iV r = 2.9 ± 0.7 x 10~ 6 ph cm" 2 s" 1 keV _i ) 
observations are also reported on the plot. 

Fig.[2]shows that we have a good agreement be- 
tween the different instruments. The HLX-1 track 
in the Hardness- Intensity Diagram (HID) is clearly 
reminiscent of that observed in GBHBs. Between 
S\ p +2p and Safter 2p, there is a softening of the 
source at nearly constant luminosity (~ 1.2 x 10 42 
erg s _1 ). Such a track in the HID is observed in 
GBHBs (see e.g. Meyer- Hofmeister et a!. 2009). 
However, the situation appears to be different be- 
cause the disk temperature near the peak of the 
outbursts PI, P2 and P3 varies at nearly constant 
luminosity (from 0.24 ± 0.02 keV to 0.18 ± 0.02 
keV - see Table O . Note that if we consider the 
3cr errors on the temperature (i.e. A\ 2 = 9 for 
one parameter of interest), then we found that all 
data have consistent fcT-values given the limited 
statistical quality of the X-ray data used. How- 
ever, when merging the Si p & S2 P spectra (5i p +2 P ) 
and the S t i & S t 2 spectra (5 t i+t2), the kT- values 
(kT = 0.2401^20 for S lp+2p and kT = Q.20±° %$ 
for Sti+ti) are no longer consistent within the 3tr 
errors. 

3.2. Investigating the L — T relation 

Fig. |3] displays the disk temperature versus the 
unabsorbed 0.2-10 keV luminosity. We fitted our 
L — kT points using a L oc T a relation where a was 
left as a free parameter. We obtained a good fit 
withlog(L) = (2.4±0.7) log T+43.7±0.5 (1 (Ter- 
rors). Feng & Kaaret (2007) showed that the lumi- 
nosity L varies as L oc T- 31±0 - 5 in the ULX NGC 
1313 X-2. Using a larger sample of ULXs Ka- 
java & Poutanen (2008) found a L oc T~ 4 relation. 
King (2009) suggested that the L oc T~ 4 relation 
is expected in the super-Eddington regime pro- 
vided that the geometrical beaming factor varies 
as b oc mT 2 . From Fig.[3j the observed trend in the 
L — kT space does not follow such a correlation. 
Mizuno et al. (2001) found that some ULXs do 



1 We also added a MEKAL component with kT = 0.44^° aI 
keV to take a possible contamination from the galaxy into 
account (see Servillat et al. 2011). Only the powerlaw 
component was used to compute the hardness ratio. 



present a L oc T 2 relation. They interpreted such 
a correlation as a possible signature for moderate 
super-Eddington accretion (see also Fukue 2000; 
K03). However, we cannot exclude a T 4 relation- 
ship at the 3cr level (see also Servillat et al. 2011). 
In Fig. [3j we also show the best-fit obtained using 
the L oc T 4 relation. 

3.3. Evidence for possible line features 

When fitting the S\ P +2 P spectrum, we found 
some strong residuals around 0.6 keV (see the top 
left panel in Fig. |4|. There may be some lower 
significance residuals at higher energy as well. We 
separately investigated the S\ p and S2 P spectra, 
and both show some residuals around 0.6 keV in 
particular for S^p- Given that the residuals are 
close to the instrumental oxygen edge, we inves- 
tigated if the residual may have an instrumental 
origin. We extracted a background spectrum using 
different parts of the field of view and encompass- 
ing the time when the source was at the peak. The 
fit of the background spectrum does not reveal any 
residuals around the O edge. We checked the im- 
ages and found no evidence of contamination by 
hot pixels. When the observations were performed 
the source location on the CCD was not close to 
the bad columns. The application of the charge 
trap correction did not change the residuals. So, 
the residuals are unlikely to be instrumental. We 
added a Gaussian line redshifted at the galaxy dis- 
tance (z — 0.0224) with an intrinsic width equal 
to zero, because this parameter is otherwise not 
constrained. The fit was improved by A\ 2 = 12.7 
for 2 d.o.f. The line centroid is equal to 0.6181°;^ 
keV and the equivalent width is EW = 93±H eV 
(see the top right panel in Fig. [4]). The disk param- 
eters do not change with kT = 0.25±0.01 keV and 
-/V = I6I4 (see Table [5] for a comparison). In the 
rest of the paper, we did not include the Gaussian 
line when fitting the S\ P +2 P spectrum, because it 
did not change the spectral parameters of the K03 
disk model considered in Section 2] significantly. 

In order to further investigate the soft X-ray 
emission at ~ 0.6 keV, we inspected the RGS data 
available from the XMM2 observation. The RGS 
1 & 2 spectra being background-dominated below 
0.5 keV and above 1.3 keV, we only considered 
this energy range in our analysis. In the 0.5-1.3 
keV energy band a total of 1230 and 1170 counts 
were collected in the RGS1 and RGS2, respectively 
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(with 520 and 605 net counts). We first used an 
absorbed diskbb plus a power-law component to 
fit the data fixing the photon index of the power- 
law component and the absorbing column density 
to the values derived from the XMM2 data (Far- 
rell et al. 2009, Servillat et al. 2011). We found 
a temperature (kT = 0.19 ± 0.02 keV) consistent 
with that derived from the EPIC data (see Ta- 
ble [5j see also Farrell et al. 2009, Servillat et al. 
2011). Although the fit is statistically acceptable 
(C-Statistic = 160.4 for 113 bins) a first inspection 
of the RGS1 & 2 residuals (see the bottom panel 
in Fig. 2]) reveals that features are present at ~ 0.6 
keV, - 0.9 keV and - 1 keV. 

We then added to the baseline continuum three 
unresolved Gaussian emission lines redshifted at 
the galaxy distance. The improvement in the fit 
due to the addition of the ~ 0.6 keV emission line 
is AC/dof = 26.2/2 (corresponding to a detec- 
tion significance > 99.9% confidence). The emis- 
sion line is detected at an energy centroid of E = 
0.642±0.002 keV, a flux of F = (1.91 ±0.6) x 10~ 5 
photons cm~ 2 s _1 and an equivalent width of 
EW = 20 ±7 eV. The energy centroid is consistent 
with that derived from the Swift-XKT data, and 
the line equivalent width seems to decrease with 
the luminosity of the continuum. An additional 
line is detected at E = 0.880±^02 keV (AC/dof 
= 12.2/2), with a flux of F = (0.93 ± 0.6) x 10~ 5 
photons cm" 2 s" 1 (EW = 27 ± 12 eV). Finally, 
the addition of the third Gaussian line at E = 
0.98 ± O.OlkeV (EW = 31 ± 22 eV) gave a much 
lower statistical significance (AC/dof = 7.0/2). 
Even if additional data are needed to further in- 
vestigate the origin of these lines, we speculate 
that the most likely identification for the ~ 0.6 
keV and ~ 0.9 keV lines is Ovm Lya and Fexvm- 
Fexix. The third line, if any, might be associated 
with Nex Lya. 

4. Investigating the accretion disk struc- 
ture 

4.1. A simple approach 

In order to investigate whether the Lx oc T~ 2 4 
relation we found in Section l3~2l could be the result 
of a slim disk surrounding the BH, we first used 
the diskpbb model instead of the diskbb model. 
The diskpbb model is a multi-temperature black- 
body disk model where local disk temperature is 



given by T(r) oc r~ p with p, a free parameter (e.g. 
Mineshige et al. 1994, Hirano et al. 1995, Watarai 
et al. 2000). When p = 0.75, the model is equiv- 
alent to a diskbb model. For p < 0.75, radial 
advection starts to become important (e.g. Fukue 
2000). Since the spectral shape of the model could 
change significantly with the p- value, we fitted the 
XMM1 and XMM2 spectra as well as those from 
Swift and Chandra. The XMM3 spectrum was not 
considered because any evidence for a disk compo- 
nent is only marginal (see Servillat et al. 2011). 
The results are summarised in Table [6] In most 
cases, the p- values are not well constrained. So, we 
cannot tell from this model whether or not radial 
advection plays a role. 

4.2. The Kawaguchi (2003) disk model 

4--2.1. Description of the model 

Kawaguchi (2003) computed disk spectra in a 
self-consistent way, taking into account the effects 
of electron scattering (opacity and disk Comp- 
tonization) and the effects of the relativistic cor- 
rection (i.e. gravitational redshift and transverse 
Doppler shift) on the disk effective temperature 
in the inner part of the accretion disk. Doppler 
boosting being not considered implies that disk 
solutions are mostly seen face-on (i = 0°). The 
model table of the disk spectra^, from sub- to 
super-Eddington accretion rates with BH masses 
of 1 — 10 3 M Q , was successfully used for analysis 
of several ULXs (e.g. Foschini et al. 2006, Vier- 
dayanti et al. 2006, Okajima et al. 2006, Yoshida 
et al. 2010). We used here a new table including 
additional computation for higher BH masses, ex- 
tending up to 10 5 M Q . The key parameters of the 
model are the BH mass (M), the accretion rate 
(M) and the viscosity parameter (a). For sub- 
Eddington accretion rates (i.e. in the case of the 
standard disk), the emergent disk spectra are a- 
insensitive. Near and at super-Eddington rates, 
the disk spectra however become a-sensitive, be- 
cause electron scattering, which is density sensi- 
tive (and density is a-sensitive), begins to affect 
the emergent spectra quite strongly. The normal- 
isation is fixed using the source distance (d = 95 
Mpc). To perform the fits, we considered the 
model option 7 that takes into account the ef- 
fects of advection, electron scattering on opacity, 

2 http: / /heasarc. nasa. gov / xanadu /xspec/ models / slimdisk. html 
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Comptonization and relativistic effects (see the 
appendix I A. II for more details about the influence 
of the different model options on the spectral pa- 
rameters from HLX-1 data). It is important to 
keep in mind that the K03 code automatically 
computes at each radius how much advection is 
present for a given accretion rate. All the models 
are computed assuming a non-rotating BH. 

4.2.2. Results 

To avoid calibration uncertainties between in- 
struments affecting the fits (see Tsujimoto et al. 
2011), we decided to fit the spectra from XMM- 
Newton, Chandra and Swift separately. 

We first fitted the Swift-XRT spectra Sx p +2p, 
S C ro and Scm by tying together the 
BH mass between them. We proceeded in the 
same way with the viscosity parameter a, while 
the accretion rate was left as an independent pa- 
rameter between the different spectra. For the 
Scro and Scr3 spectra, we added a powerlaw 
component of which the photon index was tied 
together between the two spectra. This is be- 
cause we found consistent values when fitting in- 
dividually the two spectra. When we fixed the 
a-value to the default K03 model (a = 0.01) we 
did not obtain a good fit because some spectra 
display strong residuals. In order to reproduce 
the spectra, it was necessary to leave the viscos- 
ity parameter a free. We then obtained a good 
fit (x 2 /d.o.f. = 108/104). The best-fit parame- 
ters are given in Table [Jj The derived BH mass is 
M = l-8±o:? x 10 4 Mq. Fitting the data from the 
outburst P3 using the K03 disk model, we found 
M = 1.9+Q 2 x 10 4 M Q . These estimates are con- 
sistent and well inside the IMBH mass range. We 
note that the powerlaw component was not needed 
to obtain a good fit for the St3 spectrum when 
leaving the viscosity parameter free to vary. Leav- 
ing the A^-value free to investigate the sensitiv- 
ity of the BH mass with respect to the Nh value 
did not change the spectral parameters much and 
they are still consistent within the errors with the 
values derived with Njj — 4 x 10 20 cm~ 2 . In- 
deed, we found N H = 6.5±|;| x 10 20 cm" 2 and 
M = 2.8tol x 104 M o from the outbursts P0-P2, 
while we found N H = 4.1^1^ x 10 20 cm~ 2 and 
M = 1-9±S:§ x 104 M from the outburst P3. 

Second, we fitted the Chandra spectrum. For 



the K03 disk model, we estimated for the pileup 
model a grade morphing parameter of 0.33. Note 
that this parameter does not have a strong im- 
pact on the spectral parameters derived here. The 
derived BH mass is M = 1.9 ± 0.2 x 10 4 Mq. 
We then left the -/Vy-value free and we found 
Nh = 3.7±f 9 5 x 10 20 cm" 2 and a BH mass es- 
timate (M = 1.8±o 4 x 10 4 M Q ) consistent with 
those derived from Swift. 

Even if Farrell et al. (2009) showed that the 
XMM1 spectrum is well fitted by a steep pow- 
erlaw (r ~ 3.4) suggesting that the source was 
in the steep-powerlaw state as seen in some GB- 
HBs, we decided to fit the XMM1 spectrum as 
well as the XMM2 one with the K03 disk model. 
This is because the emergent spectrum from the 
K03 disk model can have a shape strongly differ- 
ent from that of a simple MCD spectrum (see Ap- 
pendix EQ). We tied the BH mass between the 
two spectra, but we left a and M to vary indepen- 
dently. First, we keep Nh fixed at 4 x 10 20 cm~ 2 . 
For XMM2, it was necessary to add a powerlaw 
component to obtain a good fit, while the addition 
of a powerlaw for XMM1 did not improve the fit. 
The derived BH mass is then M = 1.4 ± 0.1 x 10 4 
M ; which is not consistent with the Swift-XKT 
and Chandra estimates within the 90% errors. If 
we consider the 3<r errors, then the estimates from 
the different instruments are all consistent within 
the errors. Then, in order to investigate the sen- 
sitivity of the BH mass with respect to the Nh 
value, we left free the Nh parameter. We ob- 
tained a value of Nh = 5.5lJ;J x 10 20 cm -2 and a 
BH mass of M = 1.9±J;| x 10 4 M Q . The photon 
index of the powerlaw component in both cases 
is consistent within the 90% errors between the 
XMM-Newton and Swift-XRT spectra. The best- 
fit results are given in Table [3 

For XMM1, the accretion rate (m ~ 4.4) and 
the luminosity that is very different from that de- 
rived using a powerlaw (L ~ 1.3 x 10 42 erg cm~ 2 , a 
luminosity which was observed at the peak of the 
three outbursts) are commensurate with the val- 
ues found for the XMM2 and Scro spectra (see 
Table [7J . This however implies a large value of 
a = 0.13 ±0.07, even if all a- values are consis- 
tent within the 3<r errors. If we force the viscos- 
ity parameter to the lowest possible value of the 
K03 model (a = 0.01) keeping the mass fixed at 
M = 1.8 x 10 4 Mq and N H free, we did not obtain 
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a good fit (x 2 /dof = 47.3/24). The addition of 
a powerlaw improves the fit (x 2 /dof = 23.3/22), 
but in this case the powerlaw component domi- 
nates over the disk component. In this case, we 
found similar T and Ng values to those given in 
Farrell et al. (2009). 

5. Discussion 

5.1. Constraints on the BH mass 

To be able to derive a dynamical measurement 
of the BH mass in HLX-1 will be very challeng- 
ing given the distance of the source. So, we 
have to rely on indirect estimates such as spec- 
tral fitting of the X-ray data. Zampieri & Roberts 
(2009) showed that the BH mass estimate derived 
using spectral fitting can be highly variable de- 
pending on the disk model used (see their Table 
2). In the case of HLX-1, Davis et al. (2011) 
used their advanced relativistic disk model to fit 
spectra with different shapes. They found a BH 
mass within the IMBH range (from 3000 M Q 
to 3 x 10 5 Mq) with extreme and opposite as- 
sumptions. In this paper, we used the K03 disk 
model to fit X-ray spectra of HLX-1 when the 
source was in various spectral states. We note 
that this model includes spectra for a wide range 
of accretion rates from sub- to super-Eddington 
rates. In all cases for ULXs, the fits favoured 
super- or near-Eddington accretion onto a stel- 
lar mass BH. However, the fits to HLX-1 spec- 
tra do favour an IMBH solution with a BH mass 
of M = 1.8±„j x 10 4 Mq. This is an interest- 
ing result because we get three independent and 
nonetheless consistent mass estimates from three 
different instruments (Swift-XRT, XMM-Newton 
EPIC-pn and Chandra/ ACIS) and over different 
spectral shapes. This estimate is consistent with 
the observational lower and upper limits on the 
BH mass (Farrell et al. 2009, Servillat et al. 2011). 
So, using the K03 model which assumes a non- 
spinning BH and a face-on accretion disk, we find 
that an IMBH solution is favoured with an esti- 
mated mass of - 1.8 x 10 4 M Q . The BH accretes 
at the Eddington limit (to ~ 10) and radiates close 
to the Eddington luminosity at the outburst peak 
(L Edd = 2.3tl ° 6 x 1042 er S s_1 )- This corresponds 
to an Eddington fraction of fEdd — 1-llg'g con- 
sidering a bolometric luminosity of 2.5 x 10 42 erg 
s _1 on average (see Fig. [5]). 



5.2. Comparison with the results from 
Davis et al. (2011) 

Using their BHSPEC model, Davis et al. (2011) 
studied the dependency of the BH mass in HLX-1 
with the inclination and the BH spin. We can only 
compare our results with their work for a* = 
and i — 0°. In that configuration, they found 
a consistent mass estimate for XMM2 and Chan- 
dra with logM ~ 3.8 ±0.1 (Davis 2011, private 
communication). However, no solution could be 
found in this case for the Swift spectra, because 
the fEdd = Lx I L Edd ratio was pegged to the max- 
imum permitted value in their BHSPEC model 
(i.e. fEdd — 1, Davis 2011, private communica- 
tion). From Figures 2, 5 and 6 in Davis et al. 
(2011), the best-fit contours for fixed inclination 
in the a* versus log M plots overlap well between 
the simulations for the range of allowed param- 
eters. Thus, they found a consistent mass esti- 
mate between the three instruments for i = 0° & 
a* = 0.7 with log A/ - 4.0 for Swift & Chandra 
and logM ~ 3.9 for XMM2. The mass values 
derived by Davis et al. (2011) from the XMM2 
and Chandra spectra for a* = & i — 0° are 
smaller than the one derived using the K03 disk 
model. We note that the Eddington ratios (fEdd) 
they obtained are large (and close to 1). From 
our fitting results using the K03 model, we found 
similar high f E dd values near the outburst peak. 
We believe that for such high f E dd values the ef- 
fects of advection on the emergent spectra should 
be considered. 

Evaluating the level of systematic errors of the 
K03 model by comparison with alternative mod- 
els by exploring the whole parameter space is be- 
yond the scope of the present paper and deserves 
a dedicated paper. Nevertheless, we estimated 
the level of systematic errors of the K03 model 
by comparison with the KERRBB (Li et al. 2005) 
and BHSPEC models assuming a non-spinning BH 
(a* = 0) and a face-on disk (i = 0°) in the limit 
of low mass accretion rate (i.e. a few percent of 
the Eddington limit). In that limit, the accretion 
disk is expected to be very close to the relativis- 
tic model of Novikov & Thorne (1973). We found 
that the level of systematic in the BH mass be- 
tween the K03 and KERRBB models is less than 
45% when only considering the relativistic effects 
in both models (see Appendix IA.2|) . When com- 
paring the results between the K03 and BHSPEC 
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models taking into account the relativistic effects 
and the effects due to Comptonization and elec- 
tron scattering, we found that the level of system- 
atic errors in the BH mass is less than 24% (see 
Appendix IA.2p . In all cases, the values of the ac- 
cretion rate derived using the K03, KERRBB and 
BHSPEC models agree within the errors. So, pro- 
viding similar level of systematics errors over the 
whole parameter space, we are confident that our 
BH mass and accretion rate estimates are reliable. 

5.3. The disk structure 

Fig. [5] shows a linear correlation between the 
0.01-20 keV bolomctric disk luminosity (L) and 
the accretion rate so that L = r\ rh LEdd with 
radiation efficiency rj = 0.11 ±0.03. Watarai et al. 
(2000) found a similar linear correlation for rh < 
20. This implies that advection does not play a 
significant role, and that the slim disk solutions are 
not needed. At all times during the outbursts, the 
disk aspect ratio H/R is less than 0.2, from Fig. 5 
in K03. Note that H/R is mostly independent 
from the BH mass and the viscosity parameter. 

From Fig. 1 in K03, the disk appears to be 
radiation-pressure dominated in its inner regions 
for all the m-values found. Such disks might be 
viscously and thermally unstable, and they might 
give rise to outbursts (Taam & Lin 1984; Honma et 
al. 1991, Lasota & Pelat 1991; Xue et al. 2011). 
However, Hirose et al. (2009) showed via MHD 
simulations that such disks could be thermally sta- 
ble (but see Xue et al. 2011). Most observations 
of Galactic X-ray binaries in outburst radiating 
above 10% of the Eddington luminosity for a stel- 
lar mass BH do not show any evidence for such an 
instability except maybe for GRS 1915+105 (Bel- 
loni et al. 1997; Xue et al. 2011). 

5.4. The evolution of the accretion rate 

Fig. [6] shows the evolution of the accretion rate 
through the different outbursts. At the peak of 
the three outbursts and for a few weeks (~ 3 and 
^2 — 3 weeks for the outbursts P2 and P3, respec- 
tively), the BH accretes at the Eddington limit 
(rh ~ 10). This "plateau" is seen in the X-ray 
light-curve of the outbursts P2 and P3 (see the 
right panel in Fig. [lj. Recently, Lasota et al. 
(2011) proposed that the HLX-1 outbursts may 
be due to enhanced mass-transfer rate onto a pre- 



existing accretion disk when an asymptotic giant 
branch star orbiting along an eccentric orbit with 
a period of ~ 1 year passing at periastron is tidally 
stripped. Such a mechanism is known to pro- 
duce "plateaus" in light-curves (Bath & Pringle 
1981; Esin, Lasota & Hynes 2000). The latter 
authors considered an unstable disk that became 
stable after an enhanced mass-transfer rate phase. 
In the case of HLX-1, the situation is different 
since the disk is stable. Therefore, quasi-constant 
enhanced mass-transfer rate during a given time 
interval would naturally produce a "plateau" in 
the light-curve. After the "plateau", the accre- 
tion rate drops exponentially with a decay time of 
182 ± 30 days and 92 ± 15 days for the outbursts 
PI and P2, respectively. The difference in the de- 
cay time is because there was a small re-flare event 
at the end of the outburst PI (see Fig. []}. From 
the evolution of the accretion rate over time, we 
calculated the mass accreted during the outbursts 
PI and P2 assuming that the accretion rate is the 
same throughout the disk. If it is not the case, 
then the values derived below are lower limits. We 
found an accreted mass of ~ 1.2 x 10 -4 M Q and 
~ 8.2 x 10~ 5 M Q for the outbursts PI and P2, 
respectively. 

At the end of the outbursts PI and P2 i.e. 
when the source undergoes a transition from the 
high/soft state to the low/hard state, the count 
rate drops very quickly over a timescale less than 
a month (see Fig. [1]). As shown by Servillat et 
al. (2011), any disk component in the spectrum of 
HLX-1 in the low/hard state is marginal. Given 
the XMM3 luminosity of - 2 x 10 40 erg s" 1 de- 
rived from the powerlaw component and assuming 
a radiation efficiency in the low/hard state of 0.11 
as during the outbursts, we could compute an up- 
per limit on the accretion rate in the low/hard 
state of rh < 0.09 i.e. M < 3.6 x 10~ 6 M Q yr" 1 . 
This is nicely consistent with the results obtained 
by Esin et al. (1997, 1998) using the Advection 
Dominated Accretion Flow model. According to 
that model, this must be connected with a vary- 
ing (receding) inner disc radius (e.g. Dubus et al. 
2001). Higher statistical quality of the data when 
the source transits to the low/hard state will be 
needed to further investigate that possibility. 
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5.5. The L — T relation 

In the luminosity-temperature space, we found 
a correlation L oc T~ 2 4 . Fukue (2000) and 
Kawaguchi (2003) showed that near the Edding- 
ton limit (m = 10) and above a departure from the 
L — T 4 relation is expected, due to the change of 
the radial profile of the disk temperature and the 
m-sensitivc disk colour temperature via the effects 
of electron scattering. However below this limit, 
the disk is still expected to follow the L oc T 4 re- 
lation found in the standard disk model. Apart 
from the peak of the outbursts when the accre- 
tion rate is at the Eddington limit, the accretion 
rate is below this limit. Therefore, we should ex- 
pect to have alocT 4 correlation. The behaviour 
near the peak is probably responsible for the flat- 
ter correlation. To check that, we fitted our L — kT 
points for the XMM2, S C ri, S CR 2, S CR3 , S t3 
observations (for which rh < 10 - see Table [7|) us- 
ing the relation for a standard disk. We found: 
log(L) = 4 log T + 44.95 ±0.06 (1 a error). From 
Eq. 1 in Lasota et al. (2011) giving L as a func- 
tion of the colour disk temperature at the inner 
disk radius and the BH mass, we derived the fol- 
lowing relation: log(L) ~ 4 log T + 45.1 using 
7] ~ 0.11, M = 1.8 x 10 4 M and Risco = 3i?s 
(for a non-rotating BH). This demonstrates that 
the disk in HLX-1 follows the L « T 4 relation. 

6. Conclusion 

We investigated in detail the X-ray spectral 
properties of the best IMBH candidate, HLX-1, 
using multi-epoch data collected by three differ- 
ent X-ray instruments (Swift-XRT, X MM- Newton 
EPIC-pn and Chandra). To do so, we used the 
K03 disk model in order to constrain the BH mass, 
the accretion rate and the disk structure in a non 
biased way. Indeed, this disk model has the ad- 
vantage of taking into account the effects of ra- 
dial advection and covering a wide range of ac- 
cretion rates (from sub- to super-Eddington) as- 
suming a non-rotating BH and a face-on disk in- 
clination (i = 0°). We found that our multi-epoch 
data are consistent with sub-Eddington accretion 
onto a nearly 2 x 10 4 M© BH. At the peak, the 
X-ray luminosity is near the Eddington luminos- 
ity. The derived radiation efficiency is n ~ 11%. 
The disk solution we found for HLX-1 is differ- 
ent from those derived for other ULXs using the 



same disk model (e.g. Vierdayanti et al. 2006, 
Okajima et al. 2006, Yoshida et al. 2010). In- 
deed, for other ULXs the spectral fitting favours 
super- or near Eddington accretion onto a stellar 
mass BH. Here, the disk in HLX-1 likely under- 
goes sub- and near Eddington accretion onto an 
IMBH and appears to stay relatively thin. The 
inner regions of the accretion disk are dominated 
by radiation-pressure. We also showed that the 
source globally follows the L cx T A relation. At 
the outburst peak and for a further few days, the 
BH accretes at near the Eddington limit (to ~ 10 
i.e. ~ 4 x 10~ 4 M yr _1 ). The occurence of 
this "plateau" at the outburst peaks could pro- 
vide some evidence that the outburst mechanism 
is driven by enhanced mass transfer rate onto a 
pre-existing accretion disk as proposed by Lasota 
et al. (2011). After the "plateau", the accretion 
rate decreases exponentially until the end of the 
outburst. However, it is not clear from our data 
what is the geometry of the accretion flow when 
the source transits to the low/hard state, even if 
the presence of an optically thin advection domi- 
nated flow would avoid a dramatic variation in the 
accretion rate by more than two orders of magni- 
tude from the peak to the low/hard state. Finally, 
the recent outburst starting in August 2011 gives 
more credit to the ~ 1 year recurrence timescale 
proposed by Lasota et al. (2011). Our dedicated 
multi-wavelength (radio, optical and X-ray) obser- 
vations of this outburst will shed further light on 
the nature of this unique source as well as on the 
outburst mechanism. 
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A. Appendix material 

A.l. Impact of the K03 disk model computing options on the shape of the emergent spectra, 
the BH mass and the accretion rate 

The K03 model makes use of different options to compute the emergent spectra starting from a local BB 
and then including different physical effects (advection, electron opacity, relativistic effects and Comptoniza- 
tion). First, we investigate the impact of these different options on the shape of the emergent spectra as well 
as on the M and M parameters. We take the example of the S\ p +2p spectrum as an illustration, but what 
is described below is also valid for the other spectra. 

Fig. [7] shows the modification of the S\ p +2p spectrum when selecting different model options: 1) standard 
disk & local BB; 2) slim disk & local BB; 5) slim disk including relativistic effects on the BB; 6) slim 
disk including the effects of electron scattering on opacity and relativistic effects; 7) slim disk including the 
effects of electron scattering on opacity, Comptonization and relativistic effects. Here, the term "slim disk" 
implies that the K03 code automatically computes (at each radius) how much advection is present for a 
given accretion rate. 

First, we just changed the option parameters without running a new fit for option 1 in order to see the 
effects of the different model components. When using option 1, we got for the S\ p+ 2 P spectrum a BH mass 
around 1.8 x 10 3 Mq and m ~ 32 assuming a viscosity parameter of a — 0.01. Such high m-values mean 
super-Eddington accretion. However, in this case, radial advection is expected to start playing a significant 
role. The large derived M-value also implies that electron opacity and Comptonisation should play an 
important role (see Fig. 2 in Kawaguchi 2003 for instance). This results in more spectral boosting towards 
high energies. The spectrum in red in Fig. [7] shows the global shape of the data. In order to get such a shape 
when including all the above effects in a slim disk (option 7), it is clear that the effects of both electron 
scattering on opacity and Comptonisation have to be reduced. To do so, the accretion rate (in the unit of 
the Eddington rate) has to decrease; which in turn results in an increase of the BH mass. 

To check this, we fitted the S\ p +2p spectrum for each chosen option leaving free the model parameters. 
The evolution of the model parameters are summarised in Table [8] From option 1 to option 2, there is a 
relatively small variation of the M and rh parameters (M ~ 3.4 x 10 3 M and m ~ 38). From option 2 to 
option 5 (including the relativistic effects), we obtained a significant impact on the accretion rate (rh ~ 200) 
and a slight decrease of the BH mass (M ~ 2.3 x 10 3 M®). Given the large accretion rate (rh ~ 38) derived 
in option 2, the accretion disk extends below the last stable orbit (e.g. Watarai et al. 2000). The relativistic 
effects (gravitational redshift and transverse Doppler redshift) when applying option 5 strongly suppress the 
emission coming from these regions; which in turn results in the model underestimating the data. To have a 
good match between the data and the model, the accretion rate has to be increased. From option 5 to option 
6, the spectral parameters significantly change with an increase in the BH mass by one order of magnitude 
and a decrease of the accretion rate by more than one order of magnitude. From option 6 to option 7, the 
parameters do not change significantly (M ~ 1.8 x 10 4 M Q and rh ~ 10). 

From option 5, the accretion flow is in the super-Eddington regime. The effects of Comptonisation and 
electron scattering (option 6) are likely to be significant given the large rh value (see Fig. 2 in K03). This 
would in turn result in a stronger hardening of the model than what is seen in Fig. [7] However, in order to 
obtain a good fit given the shape of the Si p +2p spectrum, the effects of both Comptonisation and electron 
scattering has to be strongly minimised. This is achieved by significantly decreasing the accretion rate (the 
accretion flow is then in the sub-Eddington regime) and therefore by significantly increasing the BH mass. 
If we force the source distance from 95 Mpc (measured distance) to an arbitrary lower value (e.g. d = 3.5 
Mpc) , such a dramatic change in the BH mass and the accretion rate is not seen because the accretion flow 
always stays in the sub-Eddington regime (see Table [9]). 
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A.2. Comparison of the K03, KERRBB & BHSPEC results in the limit of the Novikov & 
Thorne (1973) relativistic disk model 

We compare the results from the K03 disk model with those obtained with KERRBB and BHSPEC 
assuming a non-spinning BH (a* = 0) and a face-on disk (i — 0°) in the limit of low mass accretion rate 
(i.e. less than 10% of the Eddington limit). In that limit, the accretion is expected to follow the predictions 
of the relativistic disk model from Novikov & Thorne (1973). In the K03 model, this limit corresponds to 
m < 1.6 LEdd/c 2 - For the comparisons with KERRBB, we limit ourselves by setting the hardening factor 
to 1 in the model (i.e. we did not take the effects of electron scattering and Comptonisation into account). 
We also did not consider the effects of limb-darkening and self-irradiation since they are not included in 
the K03 disk model. This configuration corresponds to option 5 in the K03 disk model (only including 
the relativistic effects). For the comparisons with BHSPEC, we used the option 7 (including the effects for 
Comptonisation and electron scattering) in the K03 disk model. We used the table of BHSPEC models 
bhspcc_mass_0.01.fits[l considering BH mass up to 300 M©. 

To perform the comparisons, we proceeded in both cases as follows: 

• we fitted the 5i p +2p spectrum (i.e. when the source luminosity is maximum) using the WABS*KERRBB 
and WABS*BHSPEC model with a* = 0, i = 0° and N H = 4 x 10 20 cm" 2 for different and arbitrary 
values of the source distance D so that the m values derived were below 10% of the Eddington limit. 

• from the best-fits, we simulated a model spectrum with an exposure time of 38 ks using the Swift-~XRT 
response files and the fakeit command in Xspec. 

• we fitted the faked spectra with the WABS*K03 model with alpha fixed to 0.01, the same source 
distance D as used for the KERRBB and BHSPEC fits and the option 5 and 7 when comparing the 
results with KERRBB and BHSPEC, respectively. 

Tables [T]and [2] summarise the results of the comparisons with KERRBB and BHSPEC, respectively. From 
Table [IJ we show that there is a systematic discrepancy (less than 45% in the range considered) in the BH 
mass between the two models, while the accretion rate values appear to be consistent within the errors. From 
Table [2j we obtained a good agreement between the BHSPEC and K03 models within the 3 a errors except 
for the BH mass estimate for a distance of 2.9 Mpc. Based on that latter case, the "systematic errors" in 
the BH mass are less than 24%. We also show in Tables Q] and [2] comparisons with m values larger than 10% 
of the Eddington limit. Again, we emphasize the good agreement between the K03-KERRBB models and 
K03-BHSPEC models. So, we conclude that the K03 disk model provides reliable estimates of accretion disk 
parameters (mass and accretion rate). 



http: / /www.cita.utoronto.ca/ ~swd/xspec.html 
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Table 1: Summary of the results obtained for both the KERRBB and K03 models. The errors are given at 
the 90% confidence level. 
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- 1 ) 


42 
1.8 


52+ y 7 

1 r+0.4 
i '°-0.3 


(2.018-i) 


2. 
2. 


Mass (M ) 
m (10 19 g s~ 


- 1 ) 


29 
0.85 


39+| 

0.75 ±0.08* 


(1-3218:8?) 


1.5 
1.5 


Mass (M ) 
fa (10 19 g s" 


- 1 ) 


22 
0.49 


32+i^ 

o ™+°- 04 

u - ou -0.08 


(i-08l8:S) 



° The effects of Comptonisation, electron scattering, limb-darkening and self- irradiation were not taking 
into account. 

A We used the option 5 in the K03 model that only includes the relativistic effects, 
t Accretion rate value derived from the K03 disk model. 
* Consistent at the 3 a level. 



Table 2: Summary of the results obtained for both the BHSPEC and K03 models. The errors are given at 
the 90% confidence level. 



Distance 


Parameter 


BHSPEC 


K03 A 


m° 


(Mpc) 








LEdd/c 2 


7. 
7. 


Mass (M Q ) 
I = L x /L E dd 


234 
0.146 


237t£J 

0.20018:84 * 


(3.2018:81) 


5. 
5. 


Mass (M©) 
I 


158 
0.110 


172+23 

o.ii6l8:SS 


(i-9l8: 2 2 ) 


2.9 
2.9 


Mass (M Q ) 
I 


85 

0.069 


105+* t 

0.07318:888 


(i.i8l8:81) 


2. 
2. 


Mass (M©) 
I 


54 

0.051 


55+« 

0-06418:88^ 


(i.02i8:82 *) 



We used the option 7 in the K03 model that includes the relativistic effects as well as the effects of 
Comptonisation and electron scattering. 

° Accretion rate value derived from the K03 disk model. 

* Consistent at the 3 a level. 

t Marginally consistent at 3 a. 

* Lower error bar pegged to the lowest permitted value of the model (fa = 1). 
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Fig. 5. — Linear correlation between the 0.01- 
20 keV unabsorbed luminosity derived using the 
Kawaguchi (2003) slim disk model and the accre- 
tion rate. The luminosity was computed using the 
following relation: L bo i — An d\ F bo [ with F bo [, 
the bolometric unabsorbed flux and = 95 Mpc. 
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Fig. 6. — Evolution of the accretion rate over time. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Fig. 1.— Swift-XRT Photon Counting 0.3-10 keV light-curve of HLX-1 up to 2011-09-05. The Swift-XRT 
light-curve was obtained using the Swift-XRT light-curve generator web interface with a binning of at least 
60 counts per bin (Evans et al. 2009). Left - The vertical clotted lines mark the times of the XMM-Newton 
and Chandra observations. Right - close-up of the X-ray light-curve around the outburst peaks PI, P2 and 
P3. The zero time corresponds to the time when we first observed a re- brightening from the source (PI = 
2009-08-16 and P2 = 2010-08-29) and the peak of the outburst P3 (2011-08-15). The square on the middle 
panel is the predicted on-axis Chandra count rate into the 0.3-10 keV Swz/Z-XRT energy band. Bottom panel: 
the thin points were obtained using a lower temporal binning for data from MJD = 55788 to MJD = 55790. 
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Fig. 2. — Swift-XKF hardness versus intensity diagram. The XMM-Newton (triangles) and Chandra (square) 
points are also reported on the diagram. The arrows indicate the evolution from the peak to the end of the 
outbursts. 
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Fig. 3. — Evolution of the unabsorbed 0.01-20 keV disk luminosity with the disk temperature. The values 
for Swift (diamonds), XMM-Newton (XMM2: triangles in orange) and Chandra (square in blue) are shown 
on the plot. The thick dotted-dashed line corresponds to the best ht obtained using the L oc T A relation, 
while the thick solid line corresponds to the best fit obtained using the L oc T a relation with a ~ 2.4. The 
thick dashed line corresponds to the L oc T~ 4 prediction in the super-Eddington regime (Fukue 2000; King 
2009). 



Table 3 

LOG of the XMM-Newton, Chandra and Swift observations. 



Obs. name 


Instrument 


Obs. ID 


Start date 


End date 


Good time 








(yyyy-mm-dd) 


(yyyy-mm-dd) 


(ks) 


XMM1 


XMM-Newton EPIC-pn 


0204540201 


2004-11-23 


2004-11-23 


22 


XMM2 


XMM-Newton EPIC-pn, 


0560180901 


2008-11-28 


2008-11-28 


51 


XMM2 


XMM-Newton RGS1, RGS2 


0560180901 


2008-11-28 


2008-11-28 


51 


XMM3 


XMM-Newton EPIC-pn 


0655510201 


2010-05-14 


2010-05-15 


62 


Chandra 


Chandra ACIS 


13122 


2010-09-06 


2010-09-07 


10 




SwiftrXRT 


00031287(001-139) 


2008-10-24 


2011-10-10 


360 
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Energy (keV) Energy (keV) 



t 1 1 1 1 1 1 1 1 1 1 1 1 n — | 1 1 1 1 1 1 1 1 r 




0.7 0.8 0.9 1 1.1 

Rest Energy (keV) 



Fig. 4.— X-ray lines present in the Swift-XRT and XMM-Newton RGS data. Top - Swift-XRT spectrum 
when the source luminosity was at peak for the first and second outbursts (Si p +2p)- Left: using only an 
absorbed diskbb model. Right: Same model with the addition of a Gaussian line redshifted at the galaxy 
redshift (z = 0.0224). The dotted lines correspond to the different components of the model. Bottom - ratio 
of the RGS 1 (red) and RGS2 (black) spectra from the XMM2 observation over the model (an absorbed 
diskbb + powerlaw model) . The vertical dashed lines correspond to the centroids of the detected emission 
lines. 
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Table 4 

Nomenclature of the Swift-XRT spectra. 



Spectrum name 


Comment 


Sip 


time when the source reached its maximum luminosity 




for the first outburst PI (a count rate of CR > 0.02 cts s _1 ) 


S2p 


same but for the second outburst P2 (CR > 0.02 cts s _1 ) 


Slp+2p 


Sum of the Si p and S^p spectra 


>^ajter Zp 


~ 7 davs after the neak of the second outburst P2 (CFl ^> 02 cts s — M 


ScRO 


0.01 <CR< 0.02 cts s- 1 


ScRl 


0.015 <CR< 0.02 cts s- 1 


ScR2 


0.01 <CR< 0.015 cts s- 1 


ScR3 


0.004 <CR< 0.01 cts s- 1 


ScR4 


0.002 <CR< 0.004 cts s" 1 


SCR5 


CR < 0.002 cts s- 1 


P3- 


Outburst in August 2011 (55788 < MJD < 55809) 


P3+ 


Outburst in August 2011 (55710 < MJD < 55844) 




55788 < MJD < 55802 


'5*2 


55803 < MJD < 55822 


St3 


Rcflare event in the P3 outburst (55826 < MJD < 55844) 




Fig. 7. — Evolution of the spectral shape of the 
Kawaguchi (2003) disk model when adding dif- 
ferent physical components to an initial standard 
disk model (red). The latter model was obtained 
fitting the Si p+ 2 P spectrum. 
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Table 5 

Summary of the spectral parameters when fitting the Swift-XKT data using the diskbb 

MODEL. 



(1) 

Name* 


(2) 
kT 


(3) 
N BB 


(4) 

r 


(5) 
N r 


(6) 
HR 


(7) 
Ltot 


(8) 


(9) 
\og(L bo i) 


(10) 
X 2 /dof 




0.24 ±0.01 


20+5 






0.06 + 0.01 


12.0 + 0.4 


12.0 + 0.4 


42.2 + 0.03 


42/36 


S%p 

S2p 


0.24 ±0.02 
0.22 ±0.02 


20±« 






0.06 + 0.01 

04+ 01 


12.0 + 0.5 
13.3 + 0.7 


12.0 + 0.5 
13.3 + 0.7 


42 19+ u - 04 

AO oq+0 05 


18/21 
24/15 


Safter 2p 

SCRO 

SCR3 

ScRA 

ScR5 


0.18 ±0.02 

0.18185? 

0.17 ±0.05 


71 +b8 
' L -33 

11+36 
il -8 


o 9+0.5 
z - z -0.6 
9 9+O.5 
z - z -0.6 
9 9+O.5 
z - z -0.6 
9 9+0.5 
z - z -0.6 


2.9+1:1 

1.8 + 1.5 
2.0 + 0.7 

3 +0 - 3 

U.O_ 2 


022 +u uub 
20+ 04 

u - zu -0.06 

0.37 + 0.15 

i 9+0.9 
x - z -0.5 

1.2 + 0.9 


13.3 + 0.8 

7 K+0.3 
' -°-0.6 
q q+0.4 
J -°-0.8 

1.3 + 0.2 

0.20 + 0.06 


13.3 + 0.8 
5 6+ - 3 

°-°-0.6 

2.018:5 


42.25 + 0.07 
41.90 + 0.04 
41 46+ - 10 


54/67 
54/67 
54/67 
54/67 
54/67 


ScRl 
ScR2 


0.18 ±0.02 
0.17 ±0.02 


4 i+2b 
^ i -19 

40l?* 


2.2f 
2.2+ 


2.9 + 1.2 
3.0 ± 0.7 


u - io -0.06 
n 99+O.O6 
u - zz -0.07 


o 7 +u.9 

y - ' -1.6 

6 7 +0 - 7 


7.81^ 

4-8l° ; 3 5 


42.03 + 0.06 
41.85 + 0.06 


45/51 
45/51 


s tl 

St2 

S t3 


0.20 ±0.02 
0.20 ±0.02 
0.18 ±0.03 


44t- 


2.2* 


9 + 1.9 
z -°-2.0 


023+ 006 

U.UZO_Q QQg 

n 1O+0.13 
U - i8 -0.08 


12.4 + 2.0 
10.5+1.0 
8.6 + 1.7 


12.4 + 2.0 
10.5+1.0 
6.8 + 1.2 


42.22 + 0.07 
42.15 + 0.05 
41.98 + 0.10 


15.2/15 
17.2/22 
9.5/7 



Columns - (1) spectrum name; (2) inner disk temperature in units of keV ; (3) the black-body normal- 
isation; (4) Photon index of the powerlaw component; (5) Normalisation of the powerlaw in units of 
10 -5 photons keV -1 cm~ 2 s _1 ; (6) hardness ratio defined as the ratio of the observed 1.5-10 keV flux 
over the observed 0.3-10 keV flux; (7) unabsorbed 0.2-10 keV total luminosity in units of 10 41 erg s" 1 ; 
(8) unabsorbed 0.2-10 keV disk luminosity in units of 10 41 erg s _1 ; (9) Logarithm of the unabsorbed 
0.01-20 keV disk luminosity in units of 10 41 erg s _1 ; (10) \ 2 value and number of degrees of freedom. 

* See Table [4] for the nomenclature of the different spectra. 

' Fixed parameter. 



21 



Table 6 

Summary of the p- value derived using the diskpbb model for the XMM-Newton, Swift-XRT 

AND Chandra DATA. 



Spectrum p 



XMM1 < 0.54 

XMM2 (Ui.r 

Chandra 1, 0.57+nn7 



Sip 


0.71 ±0.2 


S2p 


0.53-1 


Safter 2p 


u - oo -0.05 


ScRQ 


0.718-i 



t The value was obtained fixing N H to 4 x 10 20 cm 2 . 
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Table 7 

Summary of the spectral parameters when fitting the XMM-Newton, Swift and Chandra data 

USING THE KAWAGUCHI (2003) DISK MODEL. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 
X 2 /dop 


Spectrum 


N H 


M 


M 


a 


r 


N r 


Ltot 


Ldisk 



XMM1 
XMM2 
XMM1 
XMM2 



5.5 
5.5 



4 
4 

+1.1 

1.0 

+1.1 

1.0 



1.4 ±0.1 
1.4 ±0.1 

-^-0.3 

1 q+i-3 

1 - a -0.3 



4.8 ±0.3 

5.9 ±0.4 
4.4t° j 
5.4±g;| 

~ FTOT" 
-2.1 
+0.7 
2.7 



11 + ' 08 
u - 11 -0.05 

01 0+ 002 

U.U±U_Q QQQ 



1.6 



0.13 ±0.07 
0.010 



0.005 
0.000 
0.00!) 
-0.015 
0.017 
-0.013 



+0.3 
-0.4 

3+0.4 
'-0.5 



1.0 ±0.4 



1.3 



+0.4 
0.6 



4.2 
5.1 
4.9 
5.7 



0.2 
-0.6 
0.1 
0.3 
+0.2 
0.6 
0.1 
0.5 

-0.4 
0.5 
0.6 



hO.2 
"-0.6 
hO.l 
-0.2 

^• y -0.6 
1+0.1 

-0.3 

Tinr 

-0.4 
0.5 
0.6 



4.2^ 

4.3: 



4.9] 



283/215 
283/215 
277/214 
277/214 



Chandra 
Chandra 



4 

3.7±j : g 



1.9 ±0.2 

1 8 + 1 - 



10.2 
10.2 



0.025 
0.023 



10.7 
10.5 



10.7 
10.5 



+0.2 

0.1 
+0.2 
-0.1 
0.2 
0.1 
0.2 
0.1 



+0.3 
2.3 
+0.3 
-1.1 

5.5 ±0.5 
3.8 ±0.4 



u - uo -0.02 
n ni r+0.005 

u.uio_ 005 

020+ 009 
u.u^u_ 006 

0.01 fixed 



+0.4 
0.7 
+4.2 
-0.3 
+0.2 
-0.4 



38/28 
37.6/27 



<Slp+2p 
Safter 2p 
ScRO 
S C R3 



1.8 
1 

1.8 
1.8 



10.1 
10.5 



1.7 
1. 



0.6 
0.8 
7+0.6 
0.8 



1.4; 

1.6 



-1.6 
-0.9 
-1.5 
1.0 



11.0 
9. 
6. 
3. 



0.7 
,+4.2 
-0.3 
,+0.3 



11.0 
9.9 
5.1 



2.3 ±0.5 



108/104 
108/104 
108/104 
108/104 



ScRl 
ScR2 


4 
4 


1.8 fixed 
1.8 fixed 


7.1 ±0.9 
4 q+0-3 


016 +U - U2S 

U.U±0_Q QQg 

niq+0 024 
u.uiy_ 009 


1 6 +(Uj 

^-O — 1.3 

1 6 +0 - 9 


1 o+2.b 
i q+2.5 


c rc+0.4 
°-°-0.6 


6.9122 

4.418:5 


46/50 
46/50 


S n 


4 


1.8 fixed 


10.9t^ 


01 +U ii 
u - ul -o.oo 






10.8 ± 1.1 


10.8 ± 1.1 


37.7/41 


St2 


4 


1.8 fixed 


9.1 ± 1.0 


O-Oielg;™ 






8.8 ±0.6 


8.8 ±0.6 


37.7/41 


St3 


4 


1.8 fixed 


6 7+ 14 

' -1.2 


0.018t8;g? 






6.7 ±0.9 


6.7 ±0.9 


37.7/41 



Columns - (1) spectrum name; (2) absorption column in units of 10 20 cm~ 2 ; (3) BH mass in units of 
10 4 M Q ; (4) accretion rate in units of with Lsdd and c, the Eddington limit and the speed of light, 
respectively; (5) viscosity parameter; (6) Photon index of the powerlaw component; (7) Normalisation 
of the powerlaw in units of 10" 5 photons keV -1 cm~ 2 s _1 ; (8) unabsorbed 0.3-10 keV total luminosity 
in units of 10 41 erg s _1 ; (9) unabsorbed 0.3-10 keV disk luminosity in units of 10 41 erg s _1 ; (10) x 2 
value and number of degrees of freedom. 

f The Swift-XRT, XMM-Newton EPIC-pn and Chandra- ACIS data were fitted separately. For a given 
instrument, we fitted together all the available spectra or a sub-set of them. 
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Table 8 

Evolution of the parameters of the Kawaguchi (2003) slim disk model for different model 
options. The fits were performed using the Si p+2 p spectrum with N h = 4 x 10 20 CM -2 . 



Option 
(1) 


M 
(2) 


m 
(3) 


a 
(4) 


chi 2 /dof 
(5) ' 


1 

2 
5 
G 
7 


1-8181 

o c+0.3 
O.O_ 2 

o q+1.0 

z -°-0.6 

29lt 4 
18 ±2 


38.0+J-4 
200t%l 

lo.i^j 


0.01 fixed 
0-82+^ 
> 0.26 
14+ 001 
n5+ 006 

U.UO_ 02 


40.4/36 
40.1/35 
51.59/35 
67.21/35 
43.4/35 



Columns - (1) Computation option from the Kawaguchi (2003) slim disk model (see Section [A. 1|) : (2) 
BH mass in units of 10 3 M Q ; (3) Accretion rate in units of LB i d ; (4) viscosity parameter; (5) value of 
the x 2 an d number of degree of freedom. 



Table 9 

Evolution of the parameters of the Kawaguchi (2003) slim disk model for different 
model options forcing the source distance to 3.5 mpc instead of the measured distance of 

95 Mpc. 



Option 
(1) 


M 
(2) 


rh 
(3) 


a 
(4) 


chi 2 /dof 
(5) ' 


1 

2 
5 
6 
7 


79ili 

58±» 
52+J 

1 00 + 13 

izz_ 10 


-i 7+0.2 
iw -0.2 

2 9 +0 - 4 

q q+0.1 
°-°-0.1 

i 4+oi 


0.01 fixed 
0.01 fixed 
0.01 fixed 
0.01 fixed 
0.01 fixed 


40.2/36 
40.3/36 
40.8/36 
40.4/36 
42/36 



Columns - (1) Computation option from the Kawaguchi (2003) slim disk model (see Section [A. 1|) : (2) 
BH mass in units of M Q ; (3) Accretion rate in units of LB i d ; (4) viscosity parameter; (5) value of the 
X 2 and number of degree of freedom. 
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